Remarks on Nitsche's functional: 
The rotationally symmetric case 



Steffen Frohlich 



Abstract 



We investigate existence and stability of rotationally symmetric critical immersions 
of variational problems of higher order which were considered in |15| and |16| . 



1 Introduction 

Let fi C K 2 be a bounded and two-fold connected domain. We consider two-dimensional immersions 

X = X(u,v) = (xLfav^aPfav^aPfav)) G C 4+a (fi,R 3 )nC°(0,R 3 ), a G (0,1), (1.1) 
with the property 

W := \X U AX V \ > in ft (1.2) 

for the surface area element W = W(u,v). Here, the indices u and v are the partial derivatives w.r.t. to 
the variables u resp. v, while A means the usual vector product in K 3 . 

Let Ki — Ki(u, v), i — 1,2, denote the principle curvatures of the surface. Then, by 

vl s k 1 (u,v)+k 2 (u,v) . 

H(u,v) :— , K(u,v) := Ki(u,v)K2(u,v) (l-o) 

we introduce its mean curvature and Gaussian curvature. We are concerned with rotationally symmetric 
critical immersions w.r.t. the variational problem 

E[X]:= JJ (a + f3H 2 - jK)W dudv — ► extr! (1.4) 
n 

with positive constants a, j3, 7 G K. 

For a ^ 0, /3, 7 = 0, the functional £[X] is proportional to the classical area functional. Critical points 
of the accessory variational problem are minimal surfaces, in this case the catenoid. In section 2 we 
investigate the stability of the catenoid solution w.r.t. small perturbations. The methods used here are 
presented in more general context in |14| . 

The case a, 7 = and (3 ^ leads to Willmore's functional for which we present selected numerical 
results. Further, we refer the reader to the textbook [57], in particular chapter 7. 

Given two coaxial circular boundary curves Fi,r2 G IR 3 with common radius R > of distance d > 0, we 
consider rotationally symmetric critical points of encouraged by the J.C.C. Nitsche's treatises |15| 
and P2||. 

2 The catenoid 
The catenary curve 

Let I := [xe,x r ] C M, |7| := \x r - x e \ > 0. A critical point / G C 2+a (I,R), a G (0, 1), of the variational 
problem 



A[f] := 2tt / f(x)yjl + f'(x) 2 dx — ► extr! (2.1) 
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is a solution of the non-linear Euler-Lagrange differential equation 

fix) f'{x) = 1 + f'(x)\ xe(xt,Xr), (2.2) 

with boundary conditions yi = f(xi) and y r = f(x r ). Note that / and /" has no zeros. With suitable 
integration constants c\ , c-i G K we get the catenary solution 

f\x) = ci cosh \J^- + c 2 J , x E I. (2.3) 

Detailed calculations can be found e.g. in jS], section 4, and in 1 , chapter IV. 

Stability of graphs 

(a) Perturbation of the catenary curve 

We investigate the stability of the catenary curve / w.r.t. small perturbations. Let ip £ C 2+ "(/,K), 
a £ (0, 1), be given such that / + if> solves also the minimal surface equation (|2.2|l . that is 

(/ + ^)(/" + O-l-(/' + ^) 2 =0 in/. (2.4) 
We arrive at the non-linear differential equation 

r{x) ~ 1W 4,1 {x) + ^ {x) = 7 ^' {x)2 ~ 7 ^W'W- (2 - 5) 

Due to //" = 1 + f' 2 > the linear differential operator on the right hand side does not obey the 
maximum principle (see Proposition 6.2 of the appendix). Thus, we assume the existence of a positive 
stability function \ e C 2+a (/,R) such that 

^' {x) ~^w^ {x)+f ~M x{x) - {) in/ ' x>0 in/ ' (2 ' 6) 

With the product trick ip — ipx we calculate 

+ 2 - 7) <p'W + (xTx + 7") ^ {x) = Hip; x) (2 ' 7) 

with the non-linear right hand side 

% X ):=^' 2 ^ W "+^-^ 2 , mx) = 0- (2.8) 

Now, the left hand side differential operator in l|2.7|) obeys the maximum principle. Corollary 6.3 ensures 
the uniqueness of the solution of the below boundary value problem I2.20fl . 

(b) Schauder norms 

Definition 2.1. We introduce the norms 



\u\\ j := max\u(x)\, \\u\\u := \\u\\ j + max \u'(x)\, 

x£l x£l 

\uh,i ■= ||M||i,/ + max|w"(x)|, 



(2.9) 



furthermore the Holder norms 



| u ( fc )(xi) - u^(x 2 )\ 

\\u\\k+a,i ■= \\u\\k,i + max j , k = 0,1,2, (2-10) 

xi,x 2 £I \Xi — X2 I 

xijtx 2 

as well as the semi norms 

[u] kJ :=max\uW(x)\, [u] k+aJ := max l - ^ -™ for k = 0,1,2. (2.11) 

xei x 1 ,x 2 ei \X\ — X2\ 

Xl=£x2 
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(c) Contraction property of $(r/; x) 
Lemma 2.1. It holds 

||$(77i;X)-^2;x)|U,J<Cl([[r ? l[[ 2 +a,I+lte||2+a,l)l|m-'?2||2 + a,I (2.12) 

for all 771 , 772 S C 2+Q (/,K) tuii/i the constant 

ci := {2ii x r 1 iu,/ + iix'V 1 /- 1 - x"r 1 iu,/(i + i/i^ Q ) 2 }(i + m 1 -") 2 ■ (2.13) 

Proof. We calculate 

*(»7i; x) - *(»»; x) = y W 2 - ^ 2 ) - y (W - - (7^ - y) fa? - % 2 ) 

= y W + - ^) - y *7iK - v'i) - y ^'(m - m) (2.14) 

x' 2 x" \ 

Ix ~ ~T) ^ m + m ^ m ~ f]2 " > ' 

Now, we estimate as follows: 

\mm;x)-®(v2;x)\\«,i 

< \\xr x \\«,i\Wi + V2L,iH-V2\\a,i 

+ ||x/- 1 |U ) /||r? 1 ||a,/||< - Aj + ||x/- 1 ||a,i||r? 2 / |U ) /||r? 1 - mLj 

+ llx'V 1 / -1 - x'T^lccAvi + m\Ui\\m - m\\ a ,i (2.15) 

< ||X/ _1 |Ul(IM|l+a,J + |M|l+a,/)||?7l ~ m\\l+a,I 

+ \\xf~ 1 \\ a ,l\\Vl\\aj\\Vl - V2h+aJ + ||x/ _1 |U,j||?72||2+a,l|hl - m\\aj 

+ Wx'^-'r 1 x'T'lkKIMk/ + IMk/)hi - V2\\ aj ■ 

Note that 

NU,<||,||o,,+ max Jg^M uii-a 

x 1 ,x 2 £l \Xl — Xi I 
ill 4 xi 

< \\r,\\ ,i+ l/l 1 "" max max 1^(01 

ii^u (2-lo) 

< Ikk/ + III 1 '" max Mr/Ill ,[x 

xi,x 2 ei 

xi^x 2 

< (l+|/r- Q )N| 1+ o,Z. 

In the same way we have 

WvWw < (1 + |/r- Q )lb?l| 2+Q ,/ , hlU,/ < (1 + lif-^IMb+a,/ . (2.17) 

We insert these inequalities into (|2.15|) to get 

ll*(»7i;x) -*(»»; x)lla,/ 

+ ||x/ _1 ||a,i(l + iJl^^llJhlla+a./llm - V2\\2+aJ (2.18) 

+ ||x/ _1 lla,i(i + ^-"fWmh+oAm - mh +a j 

+ \\x l2 x- l r l x'T'LAi + m^nWvih+cj + \\mh +a ,i)\\m - mh + «,i ■ 

This is the statement. □ 
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For tp 6 C 2+a (I, R) we define the linear differential operator 

C[p] := <p"{x) + 2 - y) V>'(*) + (y - + j) ^ (2 - 19) 

Due to (|2.6(l it obeys the maximum principle. 

Successively we will solve the boundary value problem (let $(y>) := <I>(y>; x)) 

ctn-m (220) 

ip(x^) = ^ , (p(av) = ip r 

with the non-linear right hand side Ij2.8[l and given boundary data ip(xi) = (pe and tp(x r ) = 'fir- For this 
we start with the function ^0=0 and consider the linear problems 

C[(fi k ] = $(<fk-i) in I, 

(2.21) 

<p(xg) — ipe , (p(x r ) — ip r for k = 1, 2, . . . 

From = $(0) = and the Schauder estimate 

Wfkh+aJ < C 2 \\${ipk-i)\\aj + C 3 m&x{\ip(x e )\, \<p(x r )\} (2.22) 

(see the global C 2+Q -estimate (|6.40ll together with the C°-estimatc (|6.7|) and the a priori constants 
^2,^3 G (0, +oo) following from it) we conclude 

ll^ilb+a,/ < C 3 max.{\(p(x e )\, |</?(x r )|} =: C 3 a, (2.23) 

where a := max{|y(a;f)|, |y>(a; r )|}. Using (|2.12l) we get 

\\<Pk+l - <Pkh+a,I < C 2 [|$(Vfc) - ^(Vfc-l)Ha,/ 

< CiC 2 (||^fc|| 2 +Q,/ + ||Vfe-l||2+a,/)||Vfe - Vk-llb+a,/ ■ 

For given £ G (0, 1) we choose a sufficiently small such that 



(2.24) 



Then, by induction we prove 



/"* z' 7 z' 7 

\<finh+ a ,i<aC 3 [l+ / 2 3 ) for n= 1,2,... (2.26) 



(2.27) 



1-e 

First, it holds ||yi||2+a,/ < aC 3 . Furthermore, l|2.24|l gives 

W^h+a.I < \\<P2 ~ <Plh+a,I + \\<Plh+a,I < ClC 2 \\ipi\\l +a j + \\<Plh+ a ,I 

< aC 3 (l + aCiC 2 C 3 ). 



Thus, the statement is true for n = 1 und n = 2. Let it be proved for n — 1,2, ... , to. For k < m we 
calculate (see l|2.24|l l 

ll^fc+i - Vfelb+a,/ < 2aCiC 2 C 3 ( 1 H — 2 3 I ||<p* - <pk-i\\2+ a ,i < e\\<Pk - <Pk-ih+ a ,i ■ (2.28) 



1-e 

It follows that (see 



ll^m+llb+a,! < \\<Plh+a,I + ^2\\<fk+l - <fkh+a,I 

k—1 

< \\<pi\\2+ a j + (e™- 1 + e m ~ 2 + ... + £ + 1)||V2 -yilbw (2-29) 



< aC 3 + a CiC 2 C 3 = aC 3 1 H 

1 — s V 1 — £ 
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This proves l|2.2u[l and we arrive at 

\\<Pk+i - Vkh+aj < e\\<Pk - <Pk-ih+ a ,i for all fc= 1,2,... (2.30) 

for e £ (0, 1). By Banach's fix point theorem the sequence {ipk}k=o,i,... converges uniformly in C 2+a to 
a function ip £ C 2+Q (/,R) which solves the boundary value problem (|2.2U|) . Via tp = ipx we have found 
also a solution of 1)2. 5f) with sufficiently small boundary data. 

Theorem 2.1. Let f £ C 2+a (I,M.) be a solution of the minimal surface equation 1^2.2}) and let it be 
stable in the sense of $2.b)) . For sufficiently small e > there exists ip £ C 2+a (I, M) with boundary values 
l^/H 2 ^)! ^ £ an d \^{ x r)\ < £ such that f + ip solves also y2.2)) . 



3 Variational problems of higher order 
The first variation 

Now we address the first variation of the functional £[X] from (|1.4|l . We introduce conformal parameters 
(it, v) £ £1 with the properties 

\X U \ 2 = W = \X V \ 2 , X U -X* V =Q infi. (3.1) 

We consider the perturbation 

X(u, v) := X(u, v) + s{(p(u, v)X u (u, v) + ip(u, v)X v (u, v) + x(u, v)N(u, v)}, (u, v) £ ft, (3.2) 
where e £ (— so,+sq). Let 

Z(u, v) := tp(u, v)X u (u, v) + tp(u, v)X v (u, v) + x(u, v)N(u, v). (3.3) 
From and [22] we have 

Lemma 3.1. For the first variation w.r.t. the vector field Z — Z(u,v) there hold 
5 JJ Wdudv = -2 J J H{Z ■ N^Wdudv - J[Z,N,X']ds, 

6 JJ KWdudv = - J K[Z,N,X']ds + J |aW ^ - ^ {2) {Z • N*) \ ds. 

n dn dn 

5 JJ H 2 Wdudv = JJ {AH + 2H(H 2 - K)W}(Z ■ N*) dudv 
n n 



(3.4) 



an an 
with the settings 

a (2) - — + — (— 2n£\ - — ( ii ^ 3 " 5 ^ 

ds ds \k 2 ds J ds \ k 2 ds J 

Here, k and r are the curvature and the torsion of the boundary curves, while n n and n g mean their 
normal curvature and geodesic curvature, resp., w.r.t. the surface X — X(u,v). 

The Euler-Lagrange differential equation 

f3{AH + 2H{H 2 - K)W} -2aHW = (3.6) 
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(3.7) 



has to be coupled with the natural boundary conditions 

= J{yK - f3H 2 -a}[Z,N,X']ds + J ^A {2) - (Z • N*) ds 

an an 

+ /(, 3 H- 7 A<«)^>< i , 
an 

Furthermore, we add the mean curvature differential system 

AX = 2H(X)X U A X v . (3.8) 
In the case of fixed boundary conditions there hold x\an = and [Z A X')\qu — 0, that is 

0H - 7A (1) = (3.9) 
(for a detailed discussion we refer to \T5\ an d |14|.) 

Differential equations for graph solutions 

Let A denote the Laplace-Beltrami operator w.r.t. the metric ds 2 — \X U \ 2 du 2 + 2X U - Xl dudv+ \X V \ 2 dv 2 . 
If (3 ^ we can rewrite (|3.6|) and l|3.8J) to get 

AH + 2H(H 2 - K) - ^ H = 0, WAX = 2H(X U A X v ). (3.10) 

P 

Now, assume the surface represents a graph (x, y, z(x, y)). We calculate 

WA d f l + zj \ _d_f^Zy \^9_f^Zy \ d ( 1 + 4 

- - dx{ W Zx ) dx\W Zy ) dy\W Zx J + dy\W' 



jy3 Z y) Z XX 2Z X ZyZ X y + (1 + Z x )Zyy^. 



(3.11) 



Thus, it holds 



Analogously we have 



W 3 

(1 + z 2 y )z xx - 2z x z y z xy + (1 + z 2 x )z yy = 2 H {l + z 2 x + z 2 y )i . (3.12) 



ti/ a tt % z y z xy jj (1 + Z y )(z x Z xx + %Zgj,) 1 + Z y 

W Ail = — — — il^ — H x H — - ii^ 



ZyZ XX ~\- Z X Z X y Z X Zy{Z X Z XX ~\- ZyZ X y^j Z X Zy 

w v w 3 J ~ ~W~ v 



ZyZxy ~f~ Zx^yy jj ^x^y^x^xy H - %y%yy) jj %x&y jj 

— — ±1 x -f- TT7 -o tlx t t 7" tl X y 

w w 3 w 

2z x Z xy ^ (1 + Z^){Zx^xy + z y z yy) jt , ^ j[ ^-r rr 

VF y W/3 y aa 



(3.13) 



for the mean curvature. 



The rotationally symmetric case 

We consider rotationally symmetric surfaces: The meridian function f(x) := z(x,0), where f(x) > for 
all x G [a^,av], rotates about the x-axis. We have 

z x (x, 0) = f'(x), z y (x, 0) = 0, 

1 (3.14) 



z xx {x, 0) = f"(x), z xy {x, 0) = 0, z yy (x, 0) = - 



/(*) 



G 



as well as 

H x {x,0) = H'(x), H y (x,0) = 0, 

H xx (x, 0) = H"(x), H xy {x, 0) = 0, H yy (x, 0) = 0. 
The mean curvature equation reads as 

- , 1 + /' 2 



while from 13.13fl we conclude 

fix) 



AH 



f" = 2H(1 + /'• 



/"(*) 



Investing 



l + /'(z) 2 \f(x) l + f'(x) 2 

f"(x 



f 



H'(x) 



l + f'(x)- 



H"(x). 



K(x) 



f{x){l + f'{xf} 
along the meridian curve, from (|3.1U|) and i|3.17[l we get 



H " + f (j-TTT 



H' + 2H 3 (1 + f' 2 ) + 2^H-^(l + f' 2 )H = 0. 



/ (3 



(3.15) 

(3.16) 

(3.17) 
(3.18) 

(3.19) 



4 Successive approximation 

Perturbation of the mean curvature equation 

We start with 

H = 



f" 



1 



(4.1) 



2(i + /«)i 2/yiT7^' 

Let / G C 2+a (J, R) solve the minimal surface equation. For a perturbation -0 S C 2+a (I,M.) we consider 
the mean curvature 

r+r i 



H 2(i + (/' + v) 2 ) f 2(/ + ^) v ^TT7 7 +W 

along the meridian curve of the rotationally symmetric surface / + tj). Using 

l i v 



(4.2) 



/ + V> f P 
1 



v/i + f 2 + 2/v + VTTT 2 (i + / ,2 )5 



(4.3) 



1 



3/V 



(1 + f' 2 + 2/'t/>' + V-' 2 ) § (1 + / /2 ) f (1 + / /2 ) § 



we calculate 



H = - 
2 



1 



3/V 



(l + /'2)f (! + /«)* 



l ^ 



(i + / ,2 ) ! 



(4.4) 



1 



2 V / l+7^ 

1 



+ ■ 



for 



2/ 2 v / i + r 

1 1,/ < e sufficiently small. 



f" 1 
1 + /' 2 / 



1 



/' 



2(l + /' 2 )§ 2(l + /' 2 )f 



1 3/'' 



/ 1 + /'' 
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Because / solves the minimal surface equation 1)2.2)) . we get 

1 V f" 

H = r i/f" 1 3- i/j' + J - - il> + ... (4.5) 

Lemma 4.1. Let the perturbation ip € C 2+Q (/,M) &e sufficiently small w.r.t. the C l -norm. Let f 6 
C +a (J, K) sofee i/ie minimal surface equation. Then it holds 

r - y v + j ^ = 2 (! + / /2 ) f # + *iW ( 4 - 6 ) 

ure'i/i i/ie mean curvature H along f + ip an d the non-linear term ^i(ip) which collects all super-linear 
terms for ip and its first and second derivatives. Moreover, with a stability function x from \2.f>\) . the 
product ip — ipX^ 1 satisfies 

with the non-linear right hand side 

$i(<P,H):= { ^ J > H + x^Vibx)- (4.8) 

X 

Let us denote the linear differential operator on the left hand side of 1)4. 7|) with C\. Then we consider the 
boundary value problem 

£iM = $i( Pl ff) in/, 
<p{xi) = 0, tp(x r ) = 0. 

We supplement the boundary value problem 14.13)1 for H. 

Perturbation of the Euler-Lagrange equation 

Let the variation f + tp, where ||</?||ij is sufficiently small and / solves the minimal surface equation, 
satisfy ()3.19)l . that is 



(4.10) 



+ 2 ll±^H-^-(l + f' 2 + 2f^ + ^)H. 

J+<P P 

Lemma 4.2. For the mean curvature along f + tp it holds 



H" + J r H--£(l + f 2 )H = $ 2 (tp,H) (4.11) 
/ P 



with the non-linear right hand side 



$ 2 (<y5,^) - -2(1 + / )H +(./ + tp ) \ — - - + ^ 2 + — p - — ^ ^ H 



- V + svs> + - v + w: + 1 (2/v + ^, j fl + 

where . . . means terms of higher order of ip and its derivatives. 



(4.12) 
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Let C 2 denote the linear differential operator on the left hand side of l|4.11|l . We consider the boundary 
value problem 

£ 2 [H] = <Z> 2 (<p,H) inl, 

(4-13) 

H(x e ) = 0, H(x r ) = 0. 

The homogeneous boundary conditions appear naturally in the case of rotationally symmetric Willmore 
surfaces (see |23 , section 6.3). 

Note that C 2 obeys the maximum principle if 

f" oe , „,o, la , 

y-^(l + /' 2 )<0 resp.. -<-. (4.14) 

Let us work on this assumption. 

We remark that for sufficiently small perturbations [iV'Ib+a,/ < £ die residual sums in (|4.3() and (|4.12() 
converge. For e < 1 and ||-ff||2+a,/ < 1 wc find constants C 4 ,Cs £ (0, +00) such that 

||3lfo>,fl)||a,J < C 4 (||#|| 2+ a,J + |k||i + a,i), 

( 4 - 15 ) 

\\*2(<P,H)\\a,i < C±(\\H\\l +ajI + \\ V \\ 2+a>I \\H\\ 2+a!l ) 

as well as 

pi(pi,ffi)-*i(y>a,Fa)||a,i 

< C 5 {\\H! - H 2 \\ 2+aJ + (\\(Plh+a,I + \\<P 2 \\ 2 + a ,l)\\<Pl - <p 2 \\ 2 + a ,l}, 
||$2(Vl,ffl)-*2(V2,ff2)||a,J ( 4 -!6) 

< Cg(||y)l||2+a,J + HV2H2+Q,/ + ||ifl||2+a,/ + \\H 2 \\2+ a ,l) ■■■■ 

■ ■ ■ ■ (\\<fl ~ <p 2 \\ 2 +a,I + \\Hi - H 2 \\ 2+a j). 

We will solve l|4.9f) and (|4.13[) successively: Let / be stable, and let l|4.14|l be fulfilled. We start with a 
pair (<p Q , H ). First, with the right hand side <J> 2 (<A)i Ho) we get an unique solution Hi from l|4.13|) . Now, 
we solve (|4.9f) with $ 1 ((^ ,^i) to get an unique solution tp-y. In this way we proceed with [tpx,H\). 

Let l^olb+aj, H-ffolh+a,/ < e with the above e G (0, 1). The Schauder estimates 

||^l|| 3 +a,/<C6ll*2(W),flo)||a,7, l|a+a,J < C 6 ||*i(vX>, Hi) \\ a ,i (4.17) 

with a constant C 6 G (0, +00) due to (|6.40() . I|4.7|l and l|4.11|l yield (see (|4. 15|) 
||#l|| 2 +a,I < C 4 C 6 (£ 3 + e 2 ) < 2C 4 C 6 e 2 , 

(4-18) 

lbi||2+a,/ < C 4 C 6 (||-ffi|| 2+Q ,j + £ 2 ) < C 4 C 6 (1 + 2C 4 C 6 )e 2 . 

We choose 

^"fi'^k)}' (4 - 19) 

Then there follow |-ffi||2+aj < £ and ||^i||2+a,/ < £• We continue this procedure to get 

||fl*||a+a,/<e, \\<Pkh+ a j < £ for k = 0,1,2,... (4.20) 

From l|4.16|) we deduce 

\\H k+ i - H k \\ 2+a ,i < 4eC 5 C 6 (\\H k - H k -i\\ 2+a ,i + \\tpk - Vk-ib+a,i) (4.21) 
as well as (for C5 and Cq be sufficiently large) 

W^k+l — VkWl+aJ < CsCeWHk+i — Hk\\ 2 + a ,I + ZsCsCeWifik ~ cpk-i^+a.i 

< AeClCl\\H k - fl- fc _i||2+a,j + &eClCl\\ Vk ~ v k -\\\ 2+a ,i ■ 



(4.22) 
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From the last inequalities we conclude 



||Vfe+l — <Pk\\2+aJ + \\Hk+l ~ H k \\ 2+a j 

< 10eC 2 C 2 (|K - <p k -ih+ a ,i + \\H k - H k ^\\ 2+a j) 
for k — 1,2,... Additionally we assume the smallness property 

1 

e < 



20 ClCl 



(4.23) 



(4.24) 



such that the sequence (y>fc,-fffc)fe=i,2,... is contractive. Then we find solutions ip G C 2+a (I , R) and G 
C 2+Q (/,R), H^lla+a,/ < £ an d ||H||2+aj < s, of the non-linear coupled boundary value problems Ij4.9|) 
and (|4.13|1 . As a meridian curve of the rotationally symmetric immersion with mean curvature H G 
C 2+a (I, R), we have / G C 4+a (J, R) due to well-known regularity results. 

Theorem 4.1. Let f G C 4+a (I, R) fee a stable solution of the minimal surface equation for given 

boundary values fg and f r . Let f satisfy M.14H . Then we find a solution pair (f,H) of lb)) and 
as a sufficiently small perturbation in the C 4+a -norm of the minimal surface. 

5 Numerical results 
The catenoid 

Two coaxial circular boundary curves of common radius r > span a minimal surface if the distance 
h > between the rings is sufficiently small. 

We investigate the dependence between the ratio h / r and the surface area A numerically. For this, we 
increase the distance h of the rings starting from h — 0.1 by adding successively Ah = 0.1 while keeping 
the radius r = 1.5088795 fixed (all calculations were done with Ken Brakke's Surface Evolver [5]). 

In the case h / r = 1.3256 there is no two-fold connected minimal surface. The given surface area is the 
total area of the two discs of radius r. 
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2.8292 3.7667 4.7001 5.6285 6.5506 


7.4656 


8.3723 


9.2694 




Vr 


0.7290 


0.7953 


0.8616 0.9278 0.9941 1.0604 


1.1267 


1.1929 


1.2592 


1.3256 


A 


10.1558 


11.0302 


11.8810 12.7189 13.5347 14.3377 


15.0976 


15.8227 


16.5026 


14.3250 







Vr = 0.0663 



Vr = 0.1325 



Vr = 0.1988 



Vr = 0.2651 







Vr = 0.3314 



Vr = 0.3976 



V = 0.4639 



V = 0.5302 
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tyr = 0.5965 



V = 0.6627 



V = 0.7290 



Vr = 0.7953 







V, = 0.8616 



Vr = 0.9278 



V = 0.9941 



V = 1-0604 







V = 1.1267 



V = 1.1929 



V = 1.2592 



V = 1-3256 



The Willmore catenoid 



We investigate numerically critical points of Willmore's functional (11.411 with a = 0, 7 = 0, /3 = 1. The following 
tables compare the surface area A with the Willmore energy E for different distances h (r — 1). 



h 


1.0 


1.1 


1.2 


1.3 1.4 


1.5 


1.6 


1.7 


1.8 


1.9 2.0 


A 


5.98 


6.50 


6.98 7.42 7.77 


8.13 


8.49 


8.87 


9.27 


9.68 10.11 


E 


0.00 


0.00 


0.00 0.00 0.02 


0.10 


0.22 


0.38 


0.56 


0.77 0.98 


h 


2.1 


2.2 


2.3 


2.4 


2.5 


2.6 


2.7 


2.8 


2.9 


A 


10.56 


11.02 


11.51 


12.01 


12.54 


13.10 


13.66 


14.25 


14.86 


E 


1.20 


1.43 


1.67 


1.90 


2.13 


2.36 


2.58 


2.80 


3.02 


h 


3.0 


4.0 


5.0 


6.0 


7.0 


8.0 


9.0 


10.0 


A 


15.50 


23.21 


34.15 


48.08 


65.37 


85.79 


111.07 


145.22 


E 


3.23 


5.02 


6.62 


7.56 


8.25 


8.72 


9.26 


9.71 







h = 10 



h= 11 



h = 12 



h = 13 
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J* 




4> 




■ 7 


h = 14 


h= 15 


ft = 16 


ft = 17 
















ft = 18 


ft= 19 


ft = 20 


ft = 21 
















ft = 22 


ft = 23 


ft = 24 


ft = 25 


ft = 26 




ft = 27 




ft = 28 




ft = 29 



It is conjectured that such catenoid-type critical points of Willmore's functional exist for all distances h > 0. 
Furthermore, for h sufficiently large the critical immersions seem to tend to the sphere S 2 - except from a small 
neighborhood of the boundary curves where the mean curvature H vanishes. Therefore, we expect W — > 4-7T « 12.56 
for r — > oo due to the fact that 4-7T is exactly the Willmore energy for the round sphere. 







ft = 40 



ft = 60 



ft = 80 



ft= 100 



6 Appendix: Linear ordinary differential equations 
The maximum principle 

Let u £ C 2+a (/,R), a £ (0, 1), / = [a, 6], solve the boundary value problem 

C[u](x) := u"(x) + p(x)u'(x) + q(x)u(x) — f(x) in /, u(a) = rji , u(b) = 7/2 . 
where p,g€ C a (/,R) and / € C C *(/,R). Multiplication with the positive function 

p(x) := exp J p(x) dx 

gives the Sturm-Liouville boundary value problem 

C[u] — (pu'Y + qu = f in I, u(a) = 771 , u(b) = 772 , 



(6.1a) 
(6.2) 
(6.1b) 
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where q := pq and / := pf . Note that p £ C 1 (7, R). 

In |23|. §26, Satz III we find the fundamental existence theorem: 

Proposition 6.1. (Existence and uniqueness) 

The inhomogeneous problem (6.1b) has a unique solution u £ C 2+a (J,R) if and only if the homogeneous problem 

C\u] = 0, u(a) = 0, u(b) = 0, (6.3) 

has only the trivial solution u = 0. 

The Holder regularity for it comes from the assumption / £ C a (I, R). Furthermore, the uniqueness has its origin 
in the following maximum principle which can be found in |23|. Erganzung I for §26. 

Proposition 6.2. (Maximum principle) 

Let u £ C 2 (/,R) solve (6.1b) with q < (which is equivalent toq<0). 
(i) If C[u] > and u has an inner positive maximum, then u = const, 
(ii) If C[u\ < and u has an inner negative minimum, then u = const. 

Corollary 6.1. (Uniqueness) 
Let u,v £ C 2 (7,R) solve 

C[u] = /, u(a) = r;i , u(b) = r/ 2 , 

(6.4) 

£[«] = /, «(o) = Vi , -"(b) = V2 ■ 
Assume that q < 0. Then it holds u = v in I. 

Proof. Consider w := u — v with C[w] = and w(a) = 0, w(b) = 0. The maximum principle yields w = 0. □ 

For q < the homogeneous problem has only the trivial solution, that is (6.1b) is uniquely solvable. 

Alternatively, if the homogeneous problem has non-trivial solutions, the inhomogeneous problem may have no 
solutions as the following example shows f |23|. §26): 

u" (x) + u(x) = 1, u(Q) = 0, u(tt) = 1. (6.5) 

A priori estimates 

Let po,Pi £ R and p[ £ R be constants such that 

< po < p(x) < pi < +oo and |p' (a?) | < p[ < +oo for all x £ I. (6-6) 

The maximum principle from Proposition 6.2 gives the following bounds on the C°-norm. 

Corollary 6.2. (Estimate of the C° * -norm) 
Let u £ C 2 (/,R) solve (6.1b) with q<0. Then 

Hoi <3ma, X \u(x)\ + ^^e' ,{b - a) , v := . (6.7) 



Proof. 1. Consider the positive functions 



v(x) ~ ma, X \u(x)\ + (e bt " - e^) , x £ /, (6.8) 
xedi fie^ 11 

where (ipo — pi = 1 (note that fj, > 0), such that p + p/j, > 1. It follows that 

[p(e" x )'}' = filpe^]' = nip' + pn\e^ > /le^ > ^ . (6.9) 
Because qv < we conclude 

C[v] = (p^)' + qv < -MMf^ +qv < _||/|| . ( 6 .10) 

Thus, £[v~ u] < 0. Due to the maximum principle, v — u has no local negative minimum if it is not constant. 
Because (v — u)\ei > we have v > u if v — u ^ const. 
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2. Let v — —v, then C[v] > ||/||o,/- We have C[v — u] > 0. Due to the maximum principle, v — u has no local 
positive maximum if it is not constant. Because (v — u)\gi < we have v < u ii v — u ^ const. 

3. If v — u ^ const and v — u ^ const we have 

_ \ u{x) \ _ Mhi (e M _ < < | u(x) | + MM (c * _ in /. (6 .ll) 

The statement follows for this case. 

4. Let v — u = C, C £ R. With C = u(6) - u(6) = v(b) - n 2 we get 

|C| < w(6) + 1 772 1 < 2max|u(a:)|, (6.12) 

xEdl 

and the statement follows from |m| < |t>| + |C|. Analogously, for v — u = C we get |C*j < 2 max |tt(a;)|. This 
makes the proof complete. 

□ 

In the following calculations we can replace ||«||o,/ by the above estimate if q < 0. 
Lemma 6.1. (Estimate of the C 1 -norm) 

Let u € C 2+a (I,R) solve (6.1b). Then, for all real < ^ < \ it holds 

IMIm <- ^ «/ - + (> + jjfhj * ^) M« + ^ M~, ■ CM) 

Proof. 1. First, we estimate |u'(c)| with c := 2-^. Let d := fJ.^^ with /1 < \ sufficiently small. Set x\ := c — d 
and X2 := c + d. By the mean value theorem there exists x £ [xi, X2] such that 

, = \u(x,)-u{x 2 )\ = K^)-^)! < 1 ^ 
|xi — X2 2d d 

Because |c — x\ < d we get 

c 

K( C )| = + / u"(0df| < J||«||o,7+d sup |«"(y)| 



(6.15) 



2 .... #t(fo — O) r , 



2. Integrating {pu') r = f — qu from £ = c to £ = x yields 

X X 

p(x)u'(x) = J /(Odf- y ?(0«(0de + p(c)u'(c), (6.16) 

c c 

that is, 

MM < + gi||u|l °' /) + » n ?r ^ ll U H^ + m ^" 0) ||«|| 2+ a,i • (6.17) 

2po HPo(o — a) 2po 



The statement follows. 



Lemma 6.2. (Estimate of the C° -norm) 
Let u £ C 2+a (/,R) sofoe (fi.i&j. ITien 



□ 



\\u\\ a ,i < ||M||o,/ + (6-a) 1 - a N|i,/ (6.18) 



with the above bounds on \\u\\oj and llwlli,/. 
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Lemma 6.5. (Estimate of the C 1 -norm for Poisson's equation) 
Let u G C 2+a (I,R) solve fOffi . For all real < ^ < | it holds 



(I + u)(b — a) ,, _,. 2 + a(b — a) . 

HM < ^ + - a) ' 



Proof. As in 16.151 we have 



i«(c)i < ^ Hk, + ^ - Nk ' + ^ «° 

because [u]2,i = \\F\\oj. Integrating u"(x) — F{x) from £ = c to £ = x gives 

r -i „ b — a .. _., 2 .... /i(6 — a) ,. _,, 

Mi,/ < -j- ||F|| ,z + IMk/ + ||F|kz • 

This proves the statement. 

Furthermore, [u"]2,i = \\F\\o,i yields immediately 

Lemma 6.6. (Estimate of the C 2 -norm for Poisson's equation) 
Let u G C 2+a (I,R) solve KW(l . For all real < /x < ± it holds 

ii h , 2 + (l + /x)(6-a) , 2 + ^(6 -a) 

IMkz < g ll^lkx + Mb _ fl) Hlo,/ ■ 

As in fiHSt and jg^t 

we prove 

Lemma 6.7. (Estimate of the C a -norm for Poisson's equation) 
Let u G C 2+a (I,R) solve HTM) . Then it holds 

\\u\\ a ,i < \\u\\o,i + (b - a) 1_Q [u]i,J 

with the above estimates of ||u||o,i and 

Lemma 6.8. (Estimate of the C 1+a -norm for Poisson's equation) 
Let u G C 2+a (I,M.) solve HTM) . Then it holds 

\\ll\\l+a,I < + (6— a) 1_a M2,7 

with the above bounds on [|u||i,/ and M 2 >z- 

Lemma 6.9. (Estimate of the C 2+a -norm for Poisson's equation) 
Let u G C 2+a (I,M.) solve KM) . Then it holds 

||«||a+a,i = IMk/ + [u]a+a,i = IMk/ + II^IU,/ 

with the above bound on \\u\\a,i. 

Our main result is a global estimate of the C 2+a! -norm for regular solutions of (6.1b). 

Proposition 6.3. (Global estimate of the C 2+a -norm) 
Let u G C 2+a (I,R) solve (6.1b). Let 

Pi ■= \\p\\ci , qi ■= \\q\\ a ,i , Pi := ||p||i+a,/ ■ 

Then there exist real constants C\,Ci = C\, C*2(po,Pi, Si, Pi , b — a, a) G (0, +oo) such that 

hh+*,i < Ci||/|k,j+C 2 ||M]|o,j • 

Proof. Note that 

u"(x) = -44 u » ~ 44 "00 + 44 = 
p(x) p(x) p(x) 

From 16.351 it follows that 

||«||a+a,/ = IMk/ + [u]a+a,z < IMk/ + H-Flk.z 

2 + (l + /i)(b-o) 2 + fijb-a) 

< ~ 2 \\F ||o.z + \\F\\ aJ + Kb _ a) |Mkz • 
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1. We estimate ||_F|| ,/. From 16.291 we infer 



l^lkz < ^-||/||o,7 + ^||u||x,x + J-Ho,z 
Po Po Po 



< (1 +P'l c 4)||/||o./ H (PlCs + ?l)||u||o,Z + M— C6||ll||2+a,I ■ 

PO PO PO 



We set 



such that 



We define 



(6.43) 



c 7 := — (1 +£>ic 4 ), c 8 := — (p[c 5 + qi), c 9 := — c 6 , (6.44) 
Po Po Po 



\\F\\o,i < c 7 \\f\\ j + cg\\u\\ j + nc 9 \\u\\2+ a ,i ■ (6.45) 

2. Analogously we estimate 



\\F\\ aJ < Mh± + EL \\ u \\ 1+a ,j + \\ u \\ a>1 . (6.46) 
Po Po Po 

Cio := — (l +P1C4 + qi(b — a) 1 ~ a c 4 ), 
Po 

en := ^-(pi'c 5 + gi + gi(b-a) 1 ~ Q c 5 ), (6.47) 

C12 := — (pice + qi(b - a) 1 ~ a c a ). 
po 

Together with JrHHll and t?H2l we have 

ll^lki < cioWfWci + c n \\u\\o,i + »Cl2\\u\\ 2 + a> l . (6.48) 

Now, we set 

c 13 := 2+ + f & " ^ c 7 + c 10 , 

Cl4 := 2 + (! + ;,)(* -a) C8 + 2 + f 7) + Cll , (6.49) 
2 /i(o — a) 

2 + (l + ri(b~a) _ 
C15 := C9 + C12 . 

Summarizing I6.42H . 16.451 and 16.48H we arrive at the estimate 

M|2 + a,7 < Cl3 ||/||a,J + Cl4||u||o,/ + /iCl5 1 1 U \\ 2 + a ,1 ■ (6.50) 

Finally, we choose < fi < | such that 1 — ^icis > 0. It follows that 

\H 2+a ,i < ll/Uci + ||u||o,/ • (6.51) 

1 - (IC 15 1 — fIC 15 

Setting Ci := cia(l — /zcis) -1 and C2 := 014(1 — /icis) -1 proves the statement. □ 

The presented a priori estimates are not sharp. In the foregoing application we needed only a qualitative form of 
these constants. 
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